Pressure of the 0(N) Model in 1 + 1 Dimensions 
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Abstract. The 0(N) model in 1 + 1 dimensions presents some features in common with Yang- 
Mills theories: asymptotic freedom, trace anomaly, non-petrurbative generation of a mass gap. An 
analytical approach to determine the termodynamical properties of the 0(3) model is presented and 
compared to lattice results. Here the focus is on the pressure: it is shown how to derive the pressure 
in the CJT formalism at the one-loop level by making use of the auxiliary field method. Then, the 
pressure is compared to lattice results. 
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The 0(N) model, in 1 + 1 dimensions and at nonzero temperature T, is defined by the 
following generating functional: 



J (V - exp 



P 



dx 



dxJ^b 



(1) 



whereas g is the dimensionless coupling constant, jV is a normalization constant, 

and J??o is a simple free (Euclidean) Lagrangian J£q = -d^®' d^<i> , with 4> 2 = 4> r 4>. 

4> consists of Af scalar real fields, which for future convenience we denote as <3> f = 
(<J,7t\, ... , 7liv-i). The fields are constrained by the condition 4> 2 = N/g 2 which is in- 
corporated by the delta function in Eq. (1). This nonlinear constraint enforces the ther- 
modynamics of the model on an N — 1 dimensional hypersphere and induces the inter- 
actions between the fields. This model is interesting because it shares some properties 
with Yang-Mills theories in 4 dimensions, see the review paper [1] and refs. therein: 

(i) The coupling constant g is dimensionless and the the theory is renormalizable. It 
turns out that it has a negative /3 function, thus showing asymptotic freedom. 

(ii) The model is conformal invariant at the classical level, but, just as Yang-Mills 
theories in four dimensions, at the quantum level an energy scale emerges due to 
quantum corrections and a non-perturbative mass-gap emerges (trace anomaly) [2]: 
although the Lagrangian JZq in Eq. (1) describes Af massless fields, a nonzero mass 
m = /iexp(— 2n/g 2 ) is generated dynamically due to the interactions, where /! is the 
renormalization parameter. Since the mass is non-analytic in g, it would vanish in 
perturbation theory about g = 0. 



(iii) For N = 3 the model has instantons and, at finite temperature, calorons [3]. It is 
then possible to study in a simplified framework what is the role of these nonperturbative 
field configurations in thermodynamics (for the related case in YM theories see for 
instance refs. [4, 5] and refs. therein). 

The aim of this work is to present an analytical study of the nonlinear 0(N) model 
in 1 + 1 dimensions at nonzero temperature. We compute the pressure by employing the 
CJT formalism [6] and using the auxiliary field method in the one-loop approximation 
(for other thermodynamical quantities and the two-loop calculation and technical details 
see ref. [7]). For the case /V = 3 we compare the pressure with the results of a "first- 
principles" Monte Carlo lattice calculation using the integral method. At the one-loop 
level, a good agreement is found for small temperature, but the analytic result for the 
pressure is too large when T increases. When going to the two-loop approximation, a 
considerable improvement in the high T domain is obtained, although the \ow-T part is 
slightly worsened. 



ANALYTIC APPROACH AND RESULTS 

Following ref. [8] we rewrite the partition function in Eq. (1) by using the mathemati- 
cally well-defined (i.e. convergent) form of the 5-function 

§ $2 _ = lim N / $ ae i Jo J - y 2 v s 2; ^ 

V g 2 J e^0+ J 
thus obtaining Z Q i N \ = lim e ^ + JV \ $}(X@<Pexp[— JqcIt J^dxJzf], where 

Jgf = ^dp&dn<S> + U(<S>,a) with £/(<!>,«) = l -a +|a 2 . (2) 

The quantity a is an auxiliary field serving as a Lagrange multiplier. Upon shifting the 
fields o and a by their condensates, o — > + o and a — > a® + a, a bilinear mixing 
term of the type —iaoty arises. It can be eliminated by a further shift of the a field, 
a — > a — 4i(j>a/Ne. In this way non-diagonal propagators are avoided. 

Within the so-called CJT formalism [6] the standard expression for the effective 
potential is given by 

y eff = [/(0,oo)+ £ \ [[lnGr l (k)+Dr l (k^,a )G l (k)-l]+V 2 , 

— >■ Jk 
i=a, it ,a 

where U((j),(Xo) is the tree-level potential, Di(k;<j>,OCo) are the tree-level propagators, 
Gj(k) are the full propagators. V2 contains all two-particle-irreducible self-interaction 
terms, but at one-loop level one has V% = 0. Employing the stationary conditions for the 
effective potential 5V e ff/50 = <5V e ff/<5ao = SV e s/ 5Gi(k) = (/ = <7, n, a), one gets 
the propagators G,-(fc) = — k 2 +Aff = DJ l (k\§, 0Cq) with Mf = mf and the gap equations 



h = ioco(j> + ^- [ G a {k),ioco=^- f-^+t G a (k) + {N-l) [ G K {k) 
Ne Jk Ne g L Jk Jk 
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Eliminating iOq by the previous expressions one finds the following equations for the 
condensate and the masses: 



M 



= 
2 2 



mI = mI 



4<j) 2 
~Ne 



<t> 2 -^2+ [ G a (k) + (N-l) I G n (k) 
g z Jk Jk 



(4) 



In the limit £ — > + : lim e ^ +M^ = °° — > lim e ^ + JkGa(k)/£ = 0. Thus, the gap equa- 
tions read 



= Wl , M 2 a = M 2 K + 
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-(N-l) [G n (k). (5) 



To satisfy the previous equation there are two possibilities: 
(i)M2 = 0,<M0: 



f = ^ + (N-l) fG K (k) t ^ 2 -- 2 = (N-\)f 

<5 J k 5 «/ — 



dk 1 

2^ TP 
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(ii)0=O,M^O 



M 2 G =Ml=M 2 , ~2=nJ_ 
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There is no solution to the gap equations of case (i), since the integral 



/ 



[exp[(Vk 2 /T) - 1] 
Ink 



-dk 



is divergent, whereas (j) 2 —N/g 2 is finite. Therefore, contrary to the four-dimensional 
case [8], case (ii) is the right choice. This means that in two dimensions there is no 
spontaneous symmetry breaking of the global 0(N) symmetry of the nonlinear 0(N) 
model. This reflects the Mermin-Wagner-Coleman theorem, see ref. [9], which forbids 
spontaneous breakdown of a continuous symmetry in a homogeneous system in one 
spatial dimension. 

The model needs to be regularized and then renormalized. Here we just report the 
results of these steps and we refer to ref. [7] for a detailed treatment. The renormalized 
coupling constant g ren develops a dependence on the renormalization scale /i: 
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FIGURE 1. Comparison of the analytic expression for the pressure (the solid line is the one-loop result, 
the dashed line is the two-loop result) with lattice Monte-Carlo results obtained with the integral method 
for the case N — 3. The upper right line marks the Stefan-Boltzmann limit k/3. 



thus showing that the theory is asymptotically free. At nonzero T the mass of the 
excitations is M(T). For T = we obtain 

M 2 (r = 0)=m 2 = i u 2 exp(-^-) , (6) 

V <?ren / 

which shows that a mass gap is generated and that the dilatation symmetry is broken by 
quantum fluctuations (the energy scale m is dynamically generated). Then the function 
M(T) rises linearly with T when T increases: this is an expected properties for a gas 
of quasiparticles and takes place also in the deconfined phase of Yang-Mills theories, 
e.g. ref. [10] and refs. therein. 



We now turn to the pressure of the system, which is, up to a sign, identical to the 
minimum of the effective potential: p = — V e ™ n . Its renormalized form reads 

, v M(T) 2 rdkk 2 1 M(T) 2 

p(T)=N V + N f r~, N v ; 

2gf en Jo % Q) k exp{(D k /T}- 1 87T 

(7) 

Once the pressure is known one can compute the energy density p and the trace anomaly 
using the first principle of thermodynamics: p = Tdp/dT , 9 = p — p . 

The one-loop pressure of Eq. (7) is reported in Fig. 1 and is compared with a lattice 
study of this system [7] for the case N = 3. One can notice that the result is good for 
small temperatures, but it does not match the lattice data when the temperature increases. 
This result is indeed expected: the one-loop expression in Eq. (7) reduces to T 2 Nn/6 for 
large T, which corresponds to a free gas of N non-interacting particles. This result cannot 
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be correct, because the constraint expressed by the 8 function in Eq. (1) eliminates one 
degree of freedom. 

One can go beyond the one-loop results and study the pressure at the two-loop level 
using the formalism developed in [11], see also ref. [7]. A better agreement is obtained 
for large T, although the \ow-T domain is slightly worsened. Indeed, it is possible to 
show analytically that for large T the two-loop pressure approaches the correct value 
T 2 (N- l)n/6. 



CONCLUSIONS 

In this work we have studied the 0(N) model in 1+1 dimensions at nonzero T. We have 
described how a mass gap is generated for T = 0, thus demonstrating the occurrence 
of dimensional transmutation. We have then shown explicitly the analytic results for 
the pressure at the one-loop level and briefly commented on the results at the two-loop 
level. Finally, these analytic expressions for the pressure have been compared with lattice 
results, see Fig. 1. 

More advanced analytical calculations can be performed in the future and compared to 
lattice data: in fact, this system represents a good tool to test nonperturbative approaches, 
which can then be applied to the more difficult case of Yang-Mills theories in four 
dimensions. 
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